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Abstract 

Particles crossing repeatedly the surface of a shock wave can be energized by first 
order Fermi acceleration. The linear theory is successful in describing the accelera- 
tion process as long as the pressure of the accelerated particles remains negligible 
compared to the kinetic pressure of the incoming gas (the so-called test particle 
approximation). When this condition is no longer fulfilled, the shock is modified by 
the pressure of the accelerated particles in a nonlinear way, namely the spectrum 
of accelerated particles and the shock structure determine each other. In this paper 
we present the first description of the nonlinear regime of shock acceleration when 
the shock propagates in a medium where seed particles are already present. This 
case may apply for instance to supernova shocks propagating into the interstellar 
medium, where cosmic rays are in equipartition with the gas pressure. We find that 
the appearance of multiple solutions, previously found in alternative descriptions 
of the nonlinear regime, occurs also for the case of reacceleration of seed particles. 
Moreover, for parameters of concern for supernova shocks, the shock is likely to turn 
nonlinear mainly due to the presence of the pre-existing cosmic rays, rather than 
due to the acceleration of new particles from the thermal pool. We investigate here 
the onset of the nonlinear regime for the three following cases: 1) seed particles in 
equipartition with the gas pressure; 2) particles accelerated from the thermal pool; 
3) combination of 1) and 2). 
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1 Introduction 



Diffusive acceleration at newtonian shock fronts is an extensively studied phe- 
nomenon. Detailed discussions of the current status of the investigations can 
be found in some recent excellent reviews [1-4]. While much is by now well 
understood, some issues are still subjects of much debate, for the theoretical 
and phenomenological implications that they may have. The most important 
of these is the backreaction of the accelerated particles on the shock: the vio- 
lation of the test particle approximation occurs when the acceleration process 
becomes sufficiently efficient to generate pressures of the accelerated parti- 
cles which are comparable with the incoming gas kinetic pressure. Both the 
spectrum of the particles and the structure of the shock are changed by this 
phenomenon, which is therefore intrinsically nonlinear. 

At present there are three viable approaches to account for the backreaction 
of the particles upon the shock: one is based on the ever-improving numerical 
simulations [4-10] that allow a self-consistent treatment of several effects. The 
second approach is a fluid approach, and treats cosmic rays as a relativistic 
second fluid. This class of models was proposed and discussed in [11-15]. These 
models allow one to obtain the thermodynamics of the modified shocks, but 
do not provide information about the spectrum of accelerated particles. 

The third approach is analytical and may be very helpful to understand the 

physics of the nonlinear effects in a way that sometimes is difficult to achieve 
through simulations, due to their intrinsic complexity. 

In Ref. [16] a perturbative approach was adopted, in which the pressure of 
accelerated particles was treated as a small perturbation. By construction 
this method provides an answer only for weakly modified shocks. 

An alternative approach was proposed in [17-20], based on the assumption 
that the diffusion of the particles is sufficiently energy dependent that different 
parts of the fluid are affected by particles with different average energies. The 
way the calculations are carried out implies a sort of separate solution of the 
transport equation for subrelativistic and relativistic particles, so that the two 
spectra must be somehow connected at p ~ mc a posteriori. 

Recently, in [21-23], the effects of the non-linear backreaction of accelerated 
particles on the maximum achievable energy were investigated, together with 
the effects of geometry. The maximum energy of the particles accelerated in 
supernova remnants in the presence of large acceleration efficiencies was also 
studied in [24,25]. 

The need for a practical solution of the acceleration problem in the non- linear 
regime was recognized in [26], where a simple analytical broken-power-law 
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approximation of the non-linear spectra was presented. 

Recently, some promising analytical solutions of the problem of non-linear 
shock acceleration have appeared in the literature [27-29]. These solutions 
seem to avoid many of the limitations of previous approaches. 

Numerical simulations have been instrumental to identify the dramatic effects 
of the particles backreaction: they showed that even when the fraction of par- 
ticles injected from the thermal gas is relatively small, the energy channelled 
into these few particles can be an appreciable part of the kinetic energy of 
the unshocked fluid, making the test particle approach unsuitable. The most 
visible effects of the backreaction of the accelerated particles on the shock 
appear in the spectrum of the accelerated particles, which shows a peculiar 
flattening at the highest energies. The analytical approaches reproduce well 
the basic features arising from nonlinear effects in shock acceleration. 

While several calculations exist of the nonlinear effects in the shock acceler- 
ation of quasi-monochromatic particles injected at a shock surface, there is 
no description at present of how these effects appear, if they do, when the 
shock propagates in a medium where (pre) accelerated particles already exist. 
The linear theory of this phenomenon was developed by Bell [31], but has 
never been generalized to its nonlinear extension. In fact, the backreaction 
can severely affect the process of re-energization of preaccelcratcd particles: 
Bell already showed that for strong shocks the energy content of a region 
where cosmic rays were present could be easily enhanced by a factor ~ 10 at 
each shock passage, so that cquipartition could be readily reached. In these 
conditions the backreaction of the accelerated particles should be expected. 

We report here on the first analytical treatment of the shock acceleration in 
the presence of seed nonthermal particles, with the inclusion of their nonlinear 
backreaction on the shock. Our approach is a generalization of the analytical 
method introduced in [29] to describe the nonlinear shock acceleration with 
monochromatic injection of quasi-thermal particles. In fact, we present here 
also a general calculation that accounts for both thermal particles and seed 
nonthermal particles. This case may be of interest for the study of supernova 
shocks propagating through the interstellar medium (ISM) where pressure 
balance exists between gas and cosmic rays. 

Nonlinear effects in shock acceleration of thermal particles result in the ap- 
pearance of multiple solutions in certain regions of the parameter space. This 
behaviour resembles that of critical systems, with a bifurcation occurring when 
some threshold is reached. In the case of shock acceleration, it is not easy to 
find a way of discriminating among the multiple solutions when they appear. 
Neverthless, in [30], a two fiuid approach has been used to demonstrate that 
when three solutions appear, the one with intermediate efficiency for particle 
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acceleration is unstable to corrugations in the shock structure and emission of 
acustic waves. Plausibility arguments may be put forward to justify that the 
system made of the shock plus the accelerated particles may sit at the critical 
point, but the author is not aware of any real proof that this is what hap- 
pens. The physical parameters that play a role in this approach to criticality 
are the maximum momentum achievable by the particles in the acceleration 
process, the Mach number of the shock, and the injection efficiency, namely 
the fraction of thermal particles crossing the shock that are accelerated to 
nonthermal energies. The last of them, the injection efficiency, hides a crucial 
physics problem by itself, and may play an important role in establishing the 
level of shock modification. This efficiency parameter in reality is defined by 
the microphysics of the shock and should not be a free parameter of the prob- 
lem. Unfortunately, our poor knowledge of such microphysics, in particular for 
coUisionless shocks, does not allow us to estabhsh a clear and universal con- 
nection between the injection efficiency and the macroscopic shock properties. 

The paper is structured as follows: in §2 we describe the effect of non linear- 
ity on shock acceleration and our mathematical approach to describe it. In 
particular we generalize previous calculations to the case in which seed par- 
ticles exist in the region where the shock is propagating; in §3 we describe 
the gas dynamics in the presence of a non-negligible pressure of accelerated 
particles. In §4 we describe our results, with particular attention for the onset 
of the particle backreaction and for the appearance of multiple solutions. We 
conclude in 55. 



2 Non linecir shock acceleration 

In this section we solve the diffusion-convection equation for the cosmic rays 
in the most general case in which particles are injected according to some 
function Q{x,p) and the shock propagates within a region where some cosmic 
ray distribution exists, that we denote as /oo(p), spatially homogeneous before 
the shock crossing (upstream). 

For simplicity we limit ourselves to the case of one- dimensional shocks, but 
the introduction of different geometrical effects is relatively simple, and in fact 
many of our conclusions should not be affected by geometry. 

The equation that describes the diffusive transport of particles in one dimen- 
sion is 



d_ 
dx 




— u 



df{x,p) Idu df{x,p) 



+ Q{x,p)^0, 



(1) 



dx 3 dx dp 
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where we assumed stationarity {df/dt = 0). The x axis is oriented from 
upstream to downstream, as in fig. 1. The presence of pre-existing cosmic 
rays is introduced as a boundary condition at upstream infinity, by requiring 
that f{x — —oo,p) = /oo(p)- One should keep in mind that the common 
picture of a fluid whose speed is constant (ui) until it hits the shock surface is 
not appropriate for modified shocks. In fact, in this case, the pressure of the 
accelerated particles may become large enough to slow down the fluid before 
it crosses the shock surface. Therefore in general at upstream inflnity the gas 
flows at speed Uq, different from ui (fluid speed immediately upstream of the 
shock). The two quantities are approximately equal only when the accelerated 
particles do not have dynamical relevance. The injection term is taken in the 
form Q{x,p) — Qq{p)S{x). As a flrst step, we integrate eq. 1 around x — 0, 



Upstream 




Shock 



Fig. 1. Schematic view of the shock region. 

from a; = 0~ to x = 0+, denoted in flg. 1 as points "1" and "2" respectively, 
so that the following equation can be written: 









dx 


2 


dx 



J 1 



(2) 



where Ui (■U2) is the fluid speed immediately upstream (downstream) of the 
shock and /o is the particle distribution function at the shock location. By 
requiring that the distribution function downstream is independent of the 
spatial coordinate (homogeneity), we obtain = 0, so that the boundary 
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condition at the shock can be rewritten as 



dx 



J 1 



^P^{U2 - Ui) + Qo{p). 



(3) 



We can now perform the integration of eq. (1) from x = — cxd to a; = 0~ 
(point "1"), in order to take into account the boundary condition at upstream 
infinity. Using eq. (3) we obtain 

1 dfo ^ , \ , pdu 1 du Of , , 

:^P-^{U2-Ui) -uifo + uofoc + Qo{p)+ J "■^JJ^^^J d^P'd^^ 

— oo — oo 



We can now introduce a quantity Up defined as 



Up = ui- ^ J dx^f{x,p), (5) 



-oo 



whose physical meaning is instrumental to understand the nonlinear backre- 
action of particles. The function Up is the average fluid velocity experienced 
by particles with momentum p while diffusing upstream away from the shock 
surface. In other words, the effect of the average is that, instead of a con- 
stant speed Ui upstream, a particle with momentum p experiences a spatially 
variable speed, due to the pressure of the accelerated particles that is slowing 
down the fluid. Since the diffusion coefficient is in general p-dependent, par- 
ticles with different energies feel a different compression coefficient, higher at 
higher energies if, as expected, the diffusion coefficient is an increasing function 
of momentum. 

The role of can also be explained as follows: the distribution function f{x,p) 
at a distance x from the shock surface can be written as [28] 



f{x,p) = /o(p)exp 



3D{p) 



X 

j dx'u{x') 





where q{p) = — (iln/o(p)/dlnp is the local slope of fo{p) and the diffusion 
coefficient D{p) has been assumed independent of the location x. In first ap- 
proximation, we can assume that the exponential factor remains important 
when it is of order unity, namely when its argument is much less than unity. 
We can therefore introduce a distance Xp, which is the distance at which the 
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exponential equals one. This means that 



Ui 



u— 

-/ 



dx 



du 
dx 



Ui — [ui — u{x^] ~ u{xp), 



so that the speed Up can be interpreted as the fluid speed at the point Xp 
where the particles with momentum p reverse their motion in the upstream 
fluid and return to the shock. Note that the diffusion coefficient enters the 
calculation of the distance Xp but does not enter directly the calculation of Up. 
In other words, different diffusion coefficients may move the point where the 
fluid speed is Up closer to or farther from the shock surface, but do not affect 
the value of Up. This approach is similar to that introduced in [17]. 

With the introduction of Up, eq. (4) becomes 



1 dfo 



fo 



1 dUp 



+ Mo/oo + Qo{p) = 0, 



(6) 



where we used the fact that 



d f , du „ 

p-r / dx—f 

dp J dx 



dp 



p 



dfo, . , dUp 



Eq. (6) can be written in a way that resembles more the equation for shocks 
with no particles backreaction but in the presence of seed particles with dis- 
tribution fooip) and injection function Qo{p): 



p 



dp 



Up - U2 



fo 



1 dup 



- Uofoo - Qoip) 



(7) 



The solution of this equation can be written in the following implicit form: 

/o(p) = /~(p)+/r(p) = 

^ dp3[uQf^{p) + Qo{p)] I } dp' 3 



/ 

PO 



p 



Up - U2 



exp < — 



p' Up' — U2 



1 , dUp' 



(8) 



In the case of monochromatic injection with momentum p^^j at the shock 
surface, we can write 

Qo{p)-^-^^^5{p-p,^,), (9) 
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where ngas,i is the gas density immediately upstream [x — 0~) and r] parametrizes 
the fraction of the particles crossing the shock which are going to take part in 
the acceleration process. The injection term in eq. (8) becomes 



/r(p) 



PO 



dp 3Qo{p) 

p Up-U-2 



exp 



dp' 



P' Up, - U2 



1 ,dup' 



R 



sub 



1 Anp 



3 exp 

inj 



dp' 



Pinj 



P Up 



U2 



1 ,dup' 



> . 



(10) 



Here we introduced the two quantities Rgub = Ui/u2 and Rtot = uo/u2, which 
are respectively the compression factor at the gas subshock and the total com- 
pression factor between upstream infinity and downstream. The two compres- 
sion factors would be equal in the test particle approximation. For a modified 
shock, Rtot can attain values much larger than Rsub and more in general, much 
larger than 4, which is the maximum value achievable for an ordinary strong 
non-relativistic shock. The increase of the total compression factor compared 
with the prediction for an ordinary shock is responsible for the peculiar flatten- 
ing of the spectra of accelerated particles that represents a feature of nonlinear 
effects in shock acceleration. In terms of Rgub and Rtot the density immediately 
upstream is ngas,i = {po/ mp) Rtot/ Rsub- 



In eq. (10) we can introduce a dimensionless quantity U{p) — Up/uo so that 

>. (11) 



finj/ \ _ ( ^Rsub \ Wgas,! J f dp' 3RtotU{p') 

Jo yP) — 7^ 777Z\ 7 A 3 ^^P 



RtotUip) - Ij 47rp3 ^ J p' RtotU{p' 

* Pinj 



Introducing the same formalism also for the reacceleration term in eq. (8), we 
obtain the general expression 



..X ^Rtot fdPr J fdp' 3RtotU{p') 

RtotU [p) - IJ p [ I P RtotU (p') - 1 I 



^Rsuh \ mgas,i I f dp' 3RtotU{p' 



RtotU{p) - I) At^pI^ ' \ J p' RtotU{p') -If 

' Pinj 



- ^ p T :r\ • (12) 



The solution of the problem is known if the velocity field U{p) = Up/uQ is 
known. The nonlinearity of the problem reflects in the fact that U (p) is in turn 
a function of /o as it is clear from the definition of Up. In order to solve the 
problem we need to write the equations for the thermodynamics of the system 
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including the gas, the reaccelerated cosmic rays, the cosmic rays accelerated 
from the thermal pool and the shock itself. We write and solve these equations 
in the next section. 



3 The gas dynamics of modified shocks with seed particles 

The velocity, density and thermodynamic properties of the fluid can be deter- 
mined by the mass and momentum conservation equations, with the inclusion 
of the pressure of the accelerated particles and of the preexisting cosmic rays. 
We write these equations between a point far upstream {x = —oo), where 
the fluid velocity is Mq ^i-nd the density is po = mUgas^, and the point where 
the fluid upstream velocity is Up (density pp). The index p denotes quantities 
measured at the point where the fluid velocity is Up, namely at the point Xp 
that can be reached only by particles with momentum > p. 

The mass conservation implies: 

PoUo = PpUp. (13) 

Conservation of momentum reads: 

poul + Pg,o + PCR,0 = Ppul + Pg,p + PcR,p, (14) 

where Pg Q and Pg p are the gas pressures at the points x = — oo and x = Xp 
respectively, and PcR,p is the pressure in accelerated particles at the point Xp 
(we used the symbol CR to mean cosmic rays, in the sense of accelerated par- 
ticles). The mass flow in accelerated particles has reasonably been neglected. 

Our basic assumption, similar to that used in [17], is that the diffusion is 
p-dependent and more specifically that the diffusion coefficient D{p) is an 
increasing function of p. Therefore the typical distance that a particle with 
momentum p moves away from the shock is approximately Ax ~ D{p)/up, 
larger for high energy particles than for lower energy particles ^ . As a conse- 
quence, at each given point Xp only particles with momentum larger than p 
are able to affect appreciably the fluid. Strictly speaking the validity of the 
assumption depends on how strongly the diffusion coefficient depends on the 
momentum p. 

The cosmic ray pressure at Xp is the sum of two terms: one is the pressure con- 
tributed by the adiabatic compression of the cosmic rays at upstream infinity, 

^ For the cases of interest, D{p) increases with p faster than Up does, therefore Ax 
is a monotonically increasing function of p. 
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and the second is the pressure of the particles accelerated or reaccelerated at 
the shock {Pcji{p)) and able to reach the position Xp. Since only particles with 
p can reach the point x = Xp, we can write 



^CR,p = PCR,0 — + Pcr{p) 



, Pmax 

47r 



PcrM[^] +y / dppMp)h{pl (15) 



where v{p) is the velocity of particles with momentum Pmax is the maximum 
momentum achievable in the specific situation under investigation, and 7ci? 
is the adiabatic index for the accelerated particles. 

Let us consider separately the case of a strongly modified and weakly modified 
shock, in order to determine the best choice for 7ci?- In the case of strongly 
modified shocks, we will show that most energy is piled up in the region p ~ 
Pmax ^ therefore in this case we can safely adopt ^cr = 4/3, appropriate 
for a relativistic gas. For weakly modified shocks, the accelerated particles have 
an approximately power law spectrum with a slope a. It can be shown that in 
this case ^cr = Q'/S, so that the relativistic result 7^/? = 4/3 still applies for 
q; = 4 (strong shocks). For steeper spectra (a > 4) a larger adiabatic index 
should be adopted, but in those cases the solution is basically independent 
of the choice of ^cr because of the weak backreaction of the particles. For 
the purpose of carrying out our numerical calculation we will therefore always 
take 7cR = 4/3. 

The pressure of cosmic rays at upstream infinity is simply 

, Pmax 

Att C 

PcRfi^Y J "^PP'^^^P^fM, (16) 

Pm,in 



where Pmin is some minimum momentum in the spectrum of seed particles. For 
simplicity, we assume that p^m = Pinj, namely the minimum momentum of 
the seed particles coincides with the momentum at which particles are injected 
in the shock and are accelerated. 

From eq. (14) we can see that there is a maximum distance, corresponding to 
the propagation of particles with momentum Pmax such that at larger distances 
the fluid is unaffected by the accelerated particles and Up — uq. 

We will show later that for strongly modified shocks the integral in eq. (15) is 
dominated by the region p ~ Pmax- This improves even more the validity of our 
approximation PcR,p = Pcr{> p)- This also suggests that different choices for 
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the diffusion coefficient D{p) may affect tlie value of Pmax, but at fixed Pmax 
the spectra of tlie accelerated particles should not change in a significant way. 

Assuming an adiabatic compression of the gas in the upstream region, we can 
write 

P„, = P.,{ff = P,,{^J\ (IT) 



where we used mass conservation, eq. (13). The gas pressure far upstream is 
Pg,o — Poul/ (7pMo), where 7^ is the ratio of specific heats for the gas (73 = 5/3 
for an ideal gas) and Mq is the Mach number of the fluid far upstream. 

We introduce now a parameter $,cr that quantifles the relative weight of the 
cosmic ray pressure at upstream inflnity compared with the pressure of the 
gas at the same location, $^CR = PcR,o/Pg,Q- Using this parameter and the defl- 
nition of the function U {p) , the equation for momentum conservation becomes 

+ ^^ = 0. (18) 
PoUq dp 



dU 
dp 



ICRic 



CRjj-hcR+i) 



I9 



Mi 



Mi 



u- 



-(7s+l) 



Using the deflnition of Pcr and multiplying by p, this equation becomes 



p- 



dU 
dp 



^ _ JCR^jj-(jcn+l) L[/-(7g + l) 



7. Mi 



p'vip)fo{p), (19) 



where /o depends on U{p) as written in eq. (12). Eq. (19) is therefore an 
integral-differential nonlinear equation for U (p) . The solution of this equation 
also provides the spectrum of the accelerated particles. 

The last missing piece is the connection between Rguh and Rtot, the two com- 
pression factors appearing in eq. (8). The compression factor at the gas shock 
around x — can be written in terms of the Mach number Mi of the gas 
immediately upstream through the well known expression 

R _j7^ + i)M^ .... 



On the other hand, if the upstream gas evolution is adiabatic, then the Mach 
number ai x — 0~ can be written in terms of the Mach number of the fluid at 
upstream inflnity Mq as 
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so that from the expression for Rgub we obtain 



Rtot = M, 



Ig + l 



(7, + mi, - (7. - miJ 



79 + 1 



1 

79 + 1 



(21) 



Now that an expression between Rguh and Rtot has been found, eq. (19) basi- 
caUy is an equation for Rgub, with the boundary condition that U{pmax) = 1- 
Finding the value of Rsub (and the corresponding value for Rtot) such that 
U{pmax) = 1 also provides the whole function U{p) and, through eq. (8), the 
distribution function fo{p) for the particles resulting from acceleration and 
reacceleration in the nonlinear regime. If the backreaction of the accelerated 
particles is small, the test particle solution must be recovered. 



4 The nonlinecirity of acceleration and reacceleration 



In this section we investigate the shock modification due to the backreaction 
of the accelerated particles. We spht the discussion in two parts: in §4.1 we 
consider the case of pre-existing seed particles populating the region where the 
shock is propagating, and re-energized by the shock. We find the conditions for 
the onset of the nonlinear regime in which the shock gets modified by the re- 
accelerated particles. In particular we show that the critical behaviour found 
for shock acceleration and manifesting itself through the appearance of three 
sohitions, takes place also in this case. In §4.2 we explore the more complicated 
situation in which a shock accelerates a new population of particles while 
possibly reaccelerating a pre-existing population of seed particles. 



4-1 Shocks modified by re-energized seed particles 



Here we assume that a fluid is moving with speed uo in a region where the 
temperature is Tq. The Mach number is some pre-defined value Mq, which 
can be easily related to Uq and Tq. In the region of interest we assume that a 
population of seed particles is present with energy per particle already higher 
than some injection energy necessary for the particles to feel the shock as 
a discontinuity. For simplicity we assume that these seed particles have a 
spectrum which is a power law in momentum in the form: 

/oo(p) = /oo,op-^ (22) 



In this section we explore the situation in which the shock reacceleration of 
a pre-existing population of cosmic rays can modify the shock structure, but 
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there is no acceleration of particles different from those that are already present 
as seed particles. This is equivalent to turn off the injection term [r] = in Eq. 
(12)). In the next section we discuss what happens when both components are 
present. The crucial difference between the two components is that the seed 
particles are by definition already above the threshold for injection, so that 
there is no injection efficiency that instead represents such a crucial parameter 
for the case of acceleration of particles extracted from the thermal pool. 

Having in mind the case of shocks propagating in the interstellar medium 
of our Galaxy, we consider here the case in which the infiowing gas and the 
seed particles are in pressure equilibrium, namely Pg^ = Pgr- In terms of 
the parameter ^CR = Pcr/Pq^ this implies ^CR — 1- The results can then be 
easily repeated for a generic value of ^cr- 

Eq. (19) gives the quantity U{p) as a function of the momentum p for any 
choice of Rgub ^-nd Rtot, while the acceptable solutions are those with the 
right matching conditions at Pmax- For simplicity, let us assume that the Pmax 
in eq. (22) is the same maximum momentum that particles injected at the 
shock would achieve: this value only depends on the environmental conditions 
(energy losses of the particles) and/or on the geometry of the shock, which may 
allow the escape of the particles. The momentum Pmax is the same as in the 
distribution of the pre-exisiting seed particles if, for instance, the seed particles 
have been accelerated by a shock identical to the one we are considering. In 
any case, this assumption is not needed for the validity of our conclusions and 
may be easily relaxed, it simply serves to avoid parameter proliferation. Within 
this assumption, particles re-energized at the shock are simply redistributed 
in the momentum range between pinj and Pmax- 

The solution, namely the right pair of values for the compression parame- 
ters Rgy}) and Rfot [related through eq. (21)], is obtained when the solution 
corresponding to U{pmax) = 1 is selected. 

In Fig. 2, we plot U{pmax) as a function of Rtot for a shock having Mach 
number at infinity Mq = 150. The different lines are labelled by a number rep- 
resenting the logiQ of Pmax in units of mc. In other words the parameter Pmax 
changes between lO^mc and 10^ mc. The physical solutions are found by deter- 
mining the intersections of each curve with the horizontal line corresponding 
to U{pmax) = 1- The minimum momentum of the seed particles is taken as 
lO^^mc. The intersection at Rtot — 4 is close to the well known linear solution, 
namely the solution that one would obtain in the test particle approximation. 
The energy channelled into the nonthermal particles for this solution is very 
small, and the shock remains approximately unmodified. Increasing the value 
of Pmax, the solution moves toward slightly larger values of Rtot (namely the 
shock becomes more modified). For some values of Pmax, multiple solutions 
appear. In particular three regions of Pmax c^-n be identified: 
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Fig. 2. Velocity at upstream infinity in units of uq as obtained from calculations, as 
a function of the total compression factor. The physical solutions are the ones with 
U (prnax) = 1- The plot refers to the case Mq = 150 and S^cR = 1- The different lines 
are labelled by a number representing the log^^Q of Pmax iii 

units of mc. 

1) Prnax ^ Per"*- ^^^^ region (low values of Pmax) the only solution is very 
close to the one obtained in the test particle approximation. For the values of 
the parameters in Fig. 2, fsi lO^mc. 

2) p^^ < Pmax < Pc?- ill this region (intermediate values oipmax three solutions 
appear, two of which imply a strong modification of the shock, namely an 
appreciable part of the energy flowing through the shock is converted into 
energy of the accelerated particles. For the values of the parameters in Fig. 2, 
p^^) 10^ mc. 

3) Pmax > Pct'*- this region (high values of p^aa,) the solution becomes one 
again, and the shock is always strongly modified {Rtot ^ 4). 

This critical behaviour appears also when one fixes the maximum momentum 
and uses the Mach number of the shock as the order parameter. The results 
are shown in Fig. 3, where Mach numbers between 10 and 500 have been 
considered, at fixed Pmax — lO^mc. One can see that for Mach numbers below 
~ 100 there is only one solution. For Mach numbers between ~ 100 and ~ 500 
three solutions appear, one of which roughly corresponds to an unmodified 
shock. For larger Mach numbers, only this linear solution remains, and the 
shock is always only weakly modified. The same situation is plotted in Fig. 4, 
where instead of Rfot on the x-axis there is A—Rgub, and Rgub is the compression 
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Fig. 3. Velocity at upstream infinity in units of uo as obtained from calculations as 
a funtion of the total compression factor. The physical solutions arc the ones with 
U{pmax) = 1- The plot refers to the case Pmax = 10^ mc and ^CR = 1 for the Mach 
numbers indicated on the curves. 

coefficient at the gas subshock. For unmodified shocks one expects 4 — Rgub 
0. Again three regions can be identified, in the parameter Mq. 

1) For Mo < M^^^ only one solution exists and does not necessarily correspond 
to unmodified shocks. In fact one can see that for Mo = 10 < M^W ^ iqo, the 
compression coefficient at the subshock is ~ 2, and Rtot ~ 8, for the situation 
plotted in Figs. 3 and 4. Even smaller values of Mo do not give a perfectly 
unmodified shock. 

2) For M^W < Mo < M^^^) ^j^^gg solutions appear, two of which imply a 
strong modification of the shock. For the parameters used in Figs. 3 and 4, 
Mj^) ^ 500. 

3) For Mo > M^^^ only the nearly unmodified solution exists. 

In order to emphasize the fact that the multiple solutions actually correspond 
to physical solutions with very different spectral characteristics for the accel- 
erated particles we plot in Fig. 5 the spectra for Mq = 50 (for which there 
is only one solution) and Mo = 150 (for which there are three solutions). In 
both cases we used Pmax = lO^mc. The solution corresponding to Mach num- 
ber Mo = 50 has Rgub = 2.683 and Rtot = 25.588, resulting in a total pressure 
of the accelerated particles Pcr/ PqUq — 0.8946 (solid line in Fig. 5). The three 
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Fig. 5. Spectra of re-accelerated particles for the cases Mq = 50 (single solution) and 
Mq = 150 (three solutions). The maximum momentum is fixed to Pmax = lO^mc. 

dashed lines in Fig. 5 are for Mq = 150 and represent the spectra for the 
three solutions. In numbers, these solutions (from top to bottom in Fig. 5) are 
summarized as follows: 
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first solution: Rgub — 3.0775, Rtot — 55.607, Pgr/ PqUq — 0.9611 



second solution: Rgub = 3.9839, Rtot = 14.318, Pcr/ PqUq = 0.734 

third solution: R^uh = 3.999371, Rtot = 4.2548, Pgr/ PquI = 0.06 

In the cases of strong shock modification, the asymptotic spectra for {p/mc) S> 
1 are p^f{p) ~ while the hnear theory, in case of strong shocks (Mq ^ 

1), would predict p'^f{p) ~ In the region of very low energies, the spectra 
of the reaccelerated particles tend to zero, as known from the linear theory as 
well. 

4-2 Modified shocks: the general case 

Nonhnear effects in shock acceleration were first investigated in the case when 
a fraction of the thermal gas crossing the shock surface is energized to nonther- 
mal energies. In this section we wish to apply the analytical method discussed 
in §2 and §3 to the most general case in which thermal particles are accelerated 
but seed particles may already be present in the environment. 

We start our discussion with the case of acceleration of particles from the ther- 
mal distribution, in order to show that the multiple solutions already found 
in other analytical approaches [27,28] are also obtained by adopting the ap- 
proach illustrated here. In Fig. 6, we plot the U{pmax) ^ ^ function of Rtot 
for the case in which a fraction r) of the particles (as indicated in the plots) 
crossing the shock is actually accelerated to suprathermal energies. The Mach 
number is chosen to be Mq = 150 and the maximum momentum is taken 
as Pmax = lO^mc. One can easily see that Rtot ~ 4, the value for unmodi- 
fied shocks, for small values of r) (low efficiencies), while Rtot increases above 
4 for 7] > 10~^. Three solutions appear for intermediate efficiencies, while 
the solution for the accelerated particles always predicts a strongly modified 
shock for high efficiencies, rj > 3 x 10^^. An asymptotic value of Rtot ~ 60 is 
achieved for the parameters used here. This example demonstrates that the 
critical behaviour shown to appear for high values of rj in several analytical or 
semi-analytical calculations is in fact also predicted by the approach proposed 
here. Next question to ask is however whether the presence of seed particles 
can change the critical behavior of shocks. In order to answer this question we 
consider a situation in which particles are accelerated from the thermal pool 
with efficiency 77, and at the same time the shock propagates in a medium 
where the preshock pressure in seed particles equals the gas pressure. This sit- 
uation is considered to resemble that of a supernova explosion in our Galaxy, 
where cosmic rays fill the volume remaining in quasi-equipartition with the 
gas. In order to test the critical structure of the shock we calculate U{pmax) 
as a function of the compression factor Rtot between upstream infinity and 
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Fig. 6. Velocity at upstream infinity in units of uq as obtained from calculations, as 
a function of the total compression factor. Only particles injected at the shock from 
the thermal pool are considered here with efficiency rj indicated on the curves. 

downstream. Our results are plotted m Fig. 7. The lines refer to the cases 
77 = 10-^ 10-^, 5 X 10-'', 10-^, 5 X 10-^ from bottom to top, as in Fig. 6. It 
is clear from this plot that the shock may be modified by the backreaction of 
the accelerated particles even for very low values of rj, because of the presence 
of seed particles. In other words, the nonlinearity of the shock can well be 
dominated by the presence of preaccelerated particles rather than by the ac- 
celeration of particles from the thermal pool. Clearly this depends however on 
the values of the parameters (injection momentum, 77, maximum momentum 
achievable, Mach number). The injection momentum, in principle, could be 
related to rj in order to reduce the number of free parameters of the prob- 
lem. This is possible in the assumption that the particles downstream keep a 
thermal distribution. Simulations suggest that the particles injected in the ac- 
celerator are the ones with momentum a few times the thermal momentum of 
the particles downstream, so that f] can be simply calculated by integration of 
the Maxwell-Boltzmann (MB) distribution above this minimum momentum. 
Unfortunately, we do not know whether the particle distribution downstream 
is in fact Maxwell-Boltzmann-like. Moreover, we do not know exactly what is 
the threshold to impose on the injection momentum, which is a delicate issue 
because the number of particles taking part to the acceleration process would 
result from the integral of the MB distribution over its exponentially decreas- 
ing part. Therefore we preferred here to keep r] and the injection momentum 
as separate free parameters. 
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Fig. 7. Velocity at upstream infinity in units of uq as obtained from calculations, 
as a function of the total compression factor. Both particles injected at the shock 
from the thermal pool and seed particles in pressure equilibrium with the upstream 
infinity gas pressure arc considered here. The efficiencies are as in Fig. 6. 

5 Discussion and Conclusions 



We proposed a semi-analytical approach to show that the backreaction of par- 
ticles accelerated at a shock is able to affect the shock itself, in such a way that 
the shock and the accelerated particles become parts of a nonlinear system. In 
particular, for the first time we included in this kind of calculations the seed 
particles that may be present in the region where the shock is propagating 
and that can be re-energized by the shock. 

While a test-particle approach to this problem was first presented in [31], a 
nonlinear treatment was never investigated. In the pioneering work of Ref. 
[31], it was recognized that the energy of the seed particles could be enhanced 
by about one order of magnitude at each shock passage, and that after an 
infinite number of strong shocks passing through the region, the spectrum of 
the particles would tend to the asymptotic spectrum E~^^'^. Two comments 
are in order. First, the continuous increase of the cosmic ray energy due to the 
re-energization of seed particles leads unavoidably the shock to be modified 
by the nonthermal pressure, unless one starts from an uninterestingly small 
pressure of seed particles at the beginning. Second, the fact that the spectrum 
becomes flatter than E~'^, which is the result for a strong non-relativistic 
shock, implies that most of the nonthermal energy is pushed to the highest 



19 



energies, therefore the shock again can be more easily modified. Both these 
points suggest that a nonhnear treatment of shock re-acceleration is required. 

In our Galaxy, cosmic rays are observed to be in rough equipartition with the 
gas pressure and with magnetic fields, therefore supernova shocks or shocks 
generated in other environments propagate in a medium in which the seed par- 
ticles (cosmic rays) are non-negligible. In these circumstances the non-linear 
effects may be very important. We showed here that in fact for some regions 
of the parameter space, the shock is modified mainly by the backreaction of 
the seed particles rather than by the cosmic rays accelerated at the shock from 
the thermal pool. The spectra of the re- accelerated particles have also been 
calculated. 

The interesting phenomenon of the appearance of multiple shock solutions, 
already known for the case of shock acceleration, appears also for the case of 
reacceleration. This puzzling phenomenon may suggest that the shock behaves 
as a self-regulating system settling on the critical point, as proposed in [27,28]. 
On the other hand, it is possible that the multiple solutions may be the artifact 
of some of the assumptions used in the analytical approach, in particular the 
request for time independent (stationary) solutions and the fact that the role of 
the self-generated waves on the diffusion coefficient is not taken into account. 
Further investigation, in particular in the direction of a detailed comparison 
of our results with numerical simulations of shock acceleration is required in 
order to unveil the physical meaning of the multiple solutions for modified 
shocks. 

Prom the phenomenological point of view, it would be certainly worth to study 
the implications of nonlinear shock reacceleration on the nonthermal activity 
in astrophysical environments where the effect is expected to play an important 
role, in particular in the case of supernova remnants. In particular, as pointed 
out by the referee, the suggested dominance of reaccelerated ambient seed 
particles over freshly injected particles may have serious implications for the 
spectra of secondary nuclei resulting from spallation processes. 



Acknowledgments 



The author gratefully acknowledges useful discussions with M. Baring, D. 
Ellison and M. Vietri. The author is also grateful to the referee, L.O'C. Drury 
for his very useful remarks. This work was partially supported through grant 
COFIN 2002 at the Arcetri Astrophysical Observatory. 



20 



References 

[I] L.O'C. Drury, Rep. Prog. Phys. 46 (1983) 973. 

[2] R.D. Blandford, and D. Eichler, Phys. Rep. 154 (1987) 1. 

[3] E.G. Berezhko, and G.F. Krimsky, Soviet. Phys.-Uspekhi, 12 (1988) 155. 

[4] F.C. Jones, and D.C. Ellison, Space Sci. Rev. 58 (1991) 259. 

[5] A.R. Bell, MNRAS 225 (1987) 615. 

[6] D.C. Ellison, E. Mobius, and G. Paschmann, Astrophys. J. 352 (1990) 376. 
[7] D.C. Ellison, M.G. Baring, and F.C. Jones, Astrophys. J. 453 (1995) 873. 
[8] D.C. Ellison, M.G. Baring, and F.C. Jones, Astrophys. J. 473 (1996) 1029. 
[9] H. Kang, and T.W. Jones, Astrophys. J. 476 (1997) 875. 
[10] H. Kang, T.W. Jones, U.D.J. Gieseler, Astrophys. J. 579 (2002) 337. 

[II] L.O'C Drury, and H.J. Volk, Proc. lAU Symp. 94 (1980) 363. 
[12] L.O'C Drury, and H.J. Volk, Astrophys. J. 248 (1981) 344. 

[13] L.O'C Drury, W.L Axford, and D. Summers, MNRAS 198 (1982) 833. 

[14] W.I. Axford, E. Leer, and J.F. McKenzie, Astron. & Astrophys. Ill (1982) 317. 

[15] P. Duffy, L.O'C Drury, and H. Volk, Astron. & Astrophys. 291 (1994) 613. 

[16] R. Blandford, Astrophys. J. 238 (1980) 410. 

[17] D. Eichler, Astrophys. J. 277 (1984) 429. 

[18] D. Ellison, and D. Eichler, Astrophys. J. 286 (1984) 691. 

[19] D. Eichler, Astrophys. J. 294 (1985) 40. 

[20] D.C. Ellison, and D. Eichler, Phys. Rev. Lett. 55 (1985) 2735. 

[21] E.G. Berezhko, V.K. Yelshin, and L.T. Ksenofontov, Astropart. Phys. 2 (1994) 
215. 

[22] E.G. Berezhko, L.T. Ksenofontov, and V.K. Yelshin, Nucl. Phys. B (Proc. 
Suppl.) 39A (1995) 171. 

[23] E.G. Berezhko, Astropart. Phys. 5 (1996) 367. 

[24] V.S. Ptuskin and V.N. ZirakashviU, Astron. h Astrophys. 403 (2003) 1. 
[25] V.S. Ptuskin and V.N. ZirakashviU, preprint astro-ph/0306226. 
[26] E.G. Berezhko, and D.C. Ellison, Astrophys. J. 526 (1999) 385. 



21 



[27] M.A. Malkov, Astrophys. J. 485 (1997) 638. 

[28] M.A. Malkov, P.H. Diamond and H.J. Volk, Astrophys. J. Lett. 533 (2000) 171. 

[29] P. Blasi, Astropart. Phys. 16 (2002) 429. 

[30] M. Mond and L.O'C. Drury, Astron. & Astrophys. 332 (1998) 385. 
[31] A.R. Bell, MNRAS 182 (1978) 443. 



22 



